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ABSTRACT 

Solar infrared colors provide powerful constraints on the stellar effective temperature scale, but to 
this purpose they must be measured with both accuracy and precision. We achieve this requirement 
by using line-depth ratios to derive in a mo del independent w a y the infrared colors of the Sun, 
and use the latter to test the zero-point of the iCasagrande et all |2010f ) effective temperature scale, 
confirming its accuracy. Solar colors in the widely used 2MASS JHK S and WISE Wl W2 W3 WA 
systems are provided: (V - J) = 1.198, (V - H) Q = 1.484, (V - K s ) = 1.560, (J - H) Q = 0.286, 
(J-K S ) Q = 0.362, (H-K s )q = 0.076, (V-W1)q = 1.608, (V-W2) e = 1.563, (V-W3) Q = 1.552, 
[V — W4)q = 1.604. A cross check of the effective temperatures derived implementing 2MASS or 
WISE magnitudes in the infrared flux method (IRFM) confirms that the absolute calibration of the 
two systems agree within the errors, possibly suggesting a 1% offset between the two, thus validating 
extant near and mid infrared absolute calibrations. While 2MASS magnitudes are usually well suited 
to derive T c g, we find that a number of bright, solar- like stars exhibit anomalous WISE colors. In 
most cases this effect is spurious and traceable to lower quality measurements, although for a couple 
of objects (3 ± 2% of the total sample) it might be real and hints towards the presence of warm/hot 
debris disks. 

Subject headings: techniques: photometric — Sun: fundamental parameters — stars: fundamental 
parameters 



1. INTRODUCTION 

Photometric systems and filters carry information on 
various fundamental stellar properties, such as effective 
temperature (T c g), metallicity ([Fe/H]) and surface grav- 
ity (logg). Also when studying more complex systems, 
integrated magnitudes and colors of stars can be used to 
infer properties of the underlying stellar populations, by 
interpreting observations via theoretical population syn- 
thesis models. However, stars with well known physical 
parameters and colors are needed to establish how ob- 
served photometric data must be translated into physical 
quantities and placed on an absolute scale. The absolute 
calibration of photometric systems precisely deal with 
this matter and it has a long and noble history, e specially 
in us ing solar- type stars to this purpose (e.g. Uohnson 
1965t [Wamstckeij |1981t[Campins et al.l [l9"85; Rick e et al 



20081 ). Arguably, the star with best known parameters as 
well as the most important benchmark in astrophysics is 
the Sun, but for obvious reasons it can not be observed 
with the same instruments and under the same condi- 
tions applied to distant stars, thus ma king virtually im- 
possi ble to directly measure its colors (|Stebbins &: Kronl 
[1953). 

Photometry of stars with stellar properties very sim- 
ilar to the Sun provides a way to cope with this limit, 
although it is not obvious how to identify stars satisfy- 
ing such a condition in first place. Linking photometric 
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measurements to stellar parameters is in fact the goal, 
and selecting Sun-like stars based on colors would clearly 
introduce a circular argument. On the other hand, spec- 
troscopy provides an excellent way of determining T e g, 
logg and [Fe/H] in stars, and it is routinely used to this 
purpose, although it can be heavily model-dependent. 
Nevertheless, this major limit is easily overcome when 
restricting to a purely differential analysis of stars with 
spectra largely identical to a reference one. If the lat- 
ter is solar, it is thus possible to identify the stars most 
closely resembling the Sun, the so called solar- twin fl 

Over the last few years, some of us dM clendcz ' et al.l 
2006 1 : IMelendez fc Ramire3 120071: IMelendez et al.l I2006L 
2009t iRarmrez et al.ll2009fl have conducted a systematic 
search aimed to characterise and discover the best solar- 
twins in the local ~ 100 pc volume, starting from an 
initial sample of about one hundred stars in the Hippar- 
cos catalogue chosen to be broadly consistent with be- 
ing solar-like. For each candidate, high resolution, high 
signal-to-noise observations were conducted and com- 
pared to solar reference spectra (which in fact are re- 
flected Sun-light of asteroids) obtained with the same 
instrumentation and within each observing runs (at Mc- 
D onald and Las Camp anas observatories, see Section 2 
in lRamirez et aill2012f ). 

Because the procedure adopted to identify solar-twins 
does not assume any a priori T c ff, solar- twins have al- 
ready been used to set the zero-point of the effective 
temperature scale via t he infrared flux method (IRFM, 
ICasagrande et aD l2010). The effective temperature scale 
is then a basic ingredient for measuring metallicities and, 

4 According to their increasingly similarity to the Sun, stars 
can be classified as s olar-like, solar-analogs and solar- twins 
IjCavrel de S trobcl 1996|). The term "solar-twin" does not imply 
that the stars were born together with the Sun. 
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by comparison with theoretical isochrones, to derive stel- 
lar ages. Thus, the zero-point of the effective tempera- 
ture scale directly impacts basic constraints of Galactic 
chemical evolution models (e.g the metallicity distribu- 
tion function and the age-metallicity relation) as well it 
is important to correctly interpret the Sun in a Galac- 
tic context fe.g.lNordstrom et al. 2004; Casagrand e et al.l 



l20Tfl iDatson et all 12012ft . Because of its far reaching 
implicati ons, we have contin ued to investigate this topic 
(see also iHuber et al.l [2012, for a comparison between 
the angular diameters measured by interferometry with 
those obtained via the IRFM); in particular we have 
conducted dedicated observations to overcome the ma- 
jor bottleneck in linking stellar parameters to photom- 
etry, i.e. the availability of homogeneous an d high ac- 
curacy photometric data. In iMelendez et al.1 ()2010l ) we 
have presented new Stromgren uvby observations of more 
than seventy solar-analogs and derived the colors of the 
Sun in this system, which then have been used to inves- 
tigate th e zero-point of variou s metallicity scales. Sim- 
ilarly, in IRamfrez et al.l ()2012[ ) we have presented new 
UBV(RI)c photometry of 80 solar-analogs and derived 
solar colors in the widely used Johnson- Cousins system, 
obtaining a definitive value for the long debated value of 
(B - V) Q = 0.653 ±0.003. 

In this paper we finally focus on the infrared colors 
of the Sun and the tight contraints they can provide on 
the T e g scale. In fact, even though it is possible to use 
Stromgren and Johnson-Co usins colors (jMclcndc z et al.1 
120101 : IRamfrez et al.1 12012ft . infrared ones are better 
suited to this purpose, being nearly independent on blan- 
keting and surface g ravity effects for th e spectral types 
considered here (e.g. iBessell et al.lll998l ). 

In addition to this motivation, highly standardised and 
precise infrared photometry is nowadays available from 
all-sky surveys, essentially defining new standard systems 
for the years to come: 2MASS in the near-infrared, and 
the WISE satellite in the mid-infrared. Accurate solar 
colors in these two systems are thus crucial for a numbers 
of purposes. Most importantly, the reliability at which 
infrared measurements can be used to infer stellar prop- 
erties must also be assessed: while 2MASS data can be 
confidently adopted in most cases, a number of stars seem 
to exhibit anomalous WISE colors. We find that most 
of those are artifacts which can be avoided by imposing 
more stringent observational constraints, although in a 
few cases they might be real and indicate the presence of 
debris disks. 

2. SAMPLE AND PHOTOMETRIC DATA 

Our sample cons ists of the 112 stars used in 
IRamfrez et al.1 (|2012h . from whom we have adopted V 
magnitudes and stellar parameters. The latter were de- 
rived using excitation and ionization equilibrium condi- 
tions (Ramire z et al.l 120121 and references therein), and 
because of the strictly differential analysis with respect 
to the solar reference spectrum, the impact of systematic 
errors is minimised across the limited parameter space 
covered by our sample. Although a few stars in our sam- 
ple might have faint companions, the effect seems negli- 
gible for 2MASS, but we detected a few anomalous stars 
in the WISE colors. As we discuss later, these stars were 
discarded from the analysis. 

For each star we have queried the 2MASS 



J(1.25/xm) 77(1.65^) K s (2.17/zm) and the WISE 
Wl (3.4it m) W2 (4.6^m) W3 ( Uura) W4 (22/im) cat- 
alogues (jCutri et all 120031 120121 respectively) for 
photometry. Some of the bright targets have saturated 
or unreliable 2MASS magnitudes; to retain the best 
data, in a given band we consider only observations 
having photometric quality flag "A" , read flag "1" , blend 
flag "1" (i.e. one component fit to the source), contami- 
nation and confusion flag "0" (i.e. source unaffected by 
known artifacts)0. This set of flags automatically retains 
stars with photometric uncertainty ("msigcom") in a 
given band better than 0.06 mag, and for the full sample 
mean errors are 0.02 mag in J and K s and 0.03 mag in 
H band. 

Similarly, for WISE observations we restrict our analy- 
sis to measurements consistent with being point sources 
(ext = 0, meaning that no band has a reduced \ 2 > 3 
and the source is not within 5 arcsec of a 2MASS Ex- 
tended Source Catalogue entry), unaffected by known 
artifacts (ccf = 0), quality flag "A" and variability flag 
"n" or — 5 (i.e. most likely not vari ablefl The "A" 
quality flag implies a signal-to-noise ratio higher than 
10, automatically curbing large photometric uncertain- 
ties. As detailed in the WISE Explanatory Supplement, 
the Wl — 4 channels saturate at ~ 8.1, 6.7, 3.8, —0.4 mag 
respectively, although fits to the unsaturated wings of 
the PSF allow viable magnitudes to be obtained up to 
2.0, 1.5, —3.0 and —4.0 mag. Given the brightness of our 
targets, this is never a concern as the saturated pixel 
fraction is on average 0.07 in Wl and it essentially drops 
to zero in the other bands. Finally, to further decrease 
the possibility of having spurious identifications we also 
require each WISE source to have a 2MASS point source 
counterpart associated with it. All WISE sources iden- 
tified with the above constraints have a position offset 
smaller than 3 arcsec (average 1.5) with respect to the 
target coordinate. Some 20 — 30 percent of the stars do 
fall severely apart from the main locus of the color- T e g re- 
lations, especially in W2 and W4 band (no such effect is 
visible using 2MASS), although all quality records listed 
above are fulfilled (Fig. [T]); somewhat arbitrarily we ex- 
clude stars having V-W2 > 1.65 and V-WA > 1.80 and 
we shall briefly discuss this in Section SJ Altogether, this 
set of choices limits the mean (max) error to 0.03 (0.05), 
0.02 (0.02), 0.02 (0.02), 0.07 (0.11) in Wl, W2, W3, WA 
respectively. 

3. THE SOLAR INFRARED COLORS 

For solar-type stars, V— J, V—H and V— K s are known 
to display a remarkably tight correlation with T e ff, while 
being nearly independent of other parameters such as 
[Fe/H] and log g. The strong temperature sensitivity in 
solar-type stars is due to the long wavelength baseline, 
which almost brackets the region of maximum flux in 
these stars, covering part of the spectrum with similar 
continuum opacity but differing temperature sensitivity 
to the Planck function; this argument continues to hold 
also when replacing 2MASS with WISE filters (Fig. [J). 

In the literature there are various calibrations relating 
optical/near-infrared indices to effective temperatures of 

5 See http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec2_2a.html 

6 http: / /wisc2. ipac.caltech.edu/docs/release/allsky/expsup/sec2_2a. html 

7 http: / /wise2. ipac.caltech.edu/docs/release/allsky/expsup/sec6_3d. html 
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Figure 1. Color indices versus effective temperatures obtained from the IRFM. Top panels: using only stars satisfying the photometric 
criteria discussed in Section [2] The colors of the Sun are also shown with the uncertainties derived via the IRFM (Table \3\ notice that 
in the right panel T of j q is offset by a few K for clarity purposes). Bottom panels: objects having V — W2 > 1.65 or V — W4 > 1.80 
but satisfying all other photometric and quality constraints of Section [2] are highlighted in colors. Stars from the upper panels are now 
shown with empty circles. In the bottom right panel the effective temperature and color ranges have been increased to include HD 69830 
(asterisk); also shown is HD 72905 (pentagram). For both objects WISE's photometric error bars are shown (see further discussion in 
Section |4}. 



giants and dwarfs fe.g.lRidgwav et al .1119801: 1 Alonso et all 
1996HRamirez fc Melendezll2005l;ICasagrande et al.1120101; 
Bovaiian et al.ll2012aD ; for the set of filters used in this 
work, it can be easily estimated that a change of about 
0.03 mag in solar colors implies an uncertainty of about 
50 K on the zero point of the T c s scale. Therefore, to 
check the reliability of various T c g scales, accurate and 
precise colors must be derived. This is done model inde- 
pendently in Section I3.1[ while in Section 13.21 we check 
upon the zero-point of our effective temperature scale. 

For this reason, it is important that the colors un- 
der investigation are obtained directly, without resorting 
on tran sformations between different systems. ICousinsl 
(|1987al fbl) provides relations between the Stromgren 
and Johns o n- Cousins p hotometry, and more recently, 
iBilir et all (pool l20lTh have derived an extensive set 
of color transformations relating the 2MASS, WISE and 
B V(RI)g systems . Usin g those, the solar colo rs derived 
in lMelendez et all (|2010D ; [Ramirez et all(|2012f ) and here 
are usually reproduced within 0.05 mag (with better per- 



formances when transforming from the Stromgren sys- 
tem), although certain color combinations are offset by 
as much 0.10 mag, nevertheless still consistent with the 
standard deviation of the transformations reported in 
IBilir et all (I20081 [20T1 . 

3.1. Spectral Line-Depth Ratios 

Here we use the line-depth r atios (LDRs) technique 
as described in iRamirez et al.l (|2012l ) to derive the in- 
frared colors of the Sun in a model independent way. 
Briefly, this technique exploits the fact that the ratios 
of depths of spectral line pairs with very different exci- 
tation potential are excellent T e g indicators -thus corre- 
lating well with observe d colors- ess entially independent 
of [Fe/H] and \ogg (e.g. [C^av1 [l99l. For mam sequence 
stars having v sin i < 5kms _1 , as the case of our stars, 
LDRs are also weakly depended on rotational broaden- 
ing (Biazzo et al. 2007). For each set of line pairs (from 
IKovtvukh et~aTT l20031 we measured the ratios in all stars 
of our sample and linearly fitted those ratios as function 
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Table 2 

Solar colors inferred from LDR 
Measurements . 
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of the color index under consideration (Fig. [5]), after the 
exclusion of stars not satisfying the photometric quality 
requirements discussed in Section [2j From each fit the 
standard deviation of the fit minus data residual (crm) 
was also obtained. Notice that only line pairs for which 
the color vs. LDR slope was greater than 0.3 have been 
used. Slopes shallower than this imply a lower sensitiv- 
ity, leading to larger errors in the derived solar color. 
Since the slope errors are about 0.03, this criterion is 
equivalent to a 10 a cut. An example of all line pairs 
used to derive the solar (V ~ K s ) color is given in Table 
[TJ The interpolation of those fits at the solar ratio (mea- 
sured in the reflected Sun-light of asteroids with the same 
procedure used for stars) returns the color index of the 
Sun. Since we have nine reflected Sun-light observations, 
nine solar line-depth ratio values are available for each 
line pair, resulting in nine solar colors. Columns 7 and 
8 of Table Q] provide the mean and standard deviation 
(<7g S ) of those nine values. For each color index there are 
usually about one hundred pairs available, thus making 
possible to derive extremely robust colors (Table [2]), us- 
ing the weighted mean of the values obtained from each 
line-dept h ratio, where the w eight w is 1/w — er| t + of s 
(see also IRarmrez et aLll2012fl . 

3.2. Color-T c ff Relations 

The technique presented in Section 13.11 provides an el- 
egant and model independent way of determining the 
colors of the Sun. Alternatively, it is possible to per- 
form a multiple regression of the stellar parameters 
(e.g. T G g, \ogg, [Fe/H]) relevant to a given color index 
and derive that o f the Sun by solving w i th respect to its 
parameters (e.g. IHqlmberg et al.l 120061 : IMelendez et al.l 
120101 : IRarmrez et al.l 120121) . Works focusing on the T e g 
scale adopt essentially the same approach, where the 
colors of the Sun are inferred by reverting polynomial 
colo r — [Fe/H] — T e g relations derived for dwarf star s 
(e.g. IRarmrez fc Melend cz 2005; Casagrande et al.ll2006D . 
While indices in Table [5] tightly correlate with T e ff, be- 
cause of the relatively narrow parameter space covered 
by of our stars, we have verified that they do not display 
any significant dependence on logg nor [Fe/H]. In fact, 
performing a multiple linear regression with respect to 
all three parameters or T c ff only did not improve upon 
the residual nor changed within 0.001 mag the values de- 
rived. Such a simple linear relation has also the advan- 
tage of making straightforward the connection between a 
given color index and the underlying T c g scale. Depend- 
ing on the band considered, several tens of stars survive 
the quality cuts we impose on 2MASS/WISE photom- 
etry (Section [2] and Table [3j> - Using the spectroscopic 
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Figure 2. Observed colors as a function of line-depth ratio for 
the 6119.53 A (V I), 6145.02 A (Si I) pair. The solid line is a linear 
fit to the data. 
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Table 3 

Solar colors inferred from different T of j measurements: 
spectroscopic (i.e. excitation and ionization equilibrium) and the 
IRFM implemented using Tycho2+2MASS (TY2M) or 
Johnson-Cousins+2MASS (JC2M) photometry. 



color spectroscopic IRFMty2M IRFMjc2M Altars 

(V - J) 1.207 ±0.013 1.197 ±0.013 1.185 ± 0.013 87 

(V - H) 1.499 ±0.014 1.489 ±0.013 1.475 ± 0.014 87 

(V-K s ) 1.572 ±0.013 1.563 ±0.013 1.549 ± 0.013 95 

(V -Wl) 1.620 ±0.015 1.613 ±0.014 1.601 ± 0.016 52 

(V -W2) 1.576 ±0.011 1.570 ±0.011 1.561 ± 0.012 56 

(V - W3) 1.564 ±0.013 1.558 ±0.014 1.545 ± 0.016 59 

(V -WA) 1.610 ±0.021 1.601 ±0.018 1.592 ± 0.018 42 



(excitation and ionization equilibrium) temperatures de- 
termined for the full sample of solar-analogs returns solar 
indices systematically redder by ~ 0.01 mag with respect 
to those obtained via LDRs. This implies that on average 
our spectroscopic T e ff are overesti mated by about 20 K , 
in agreement with what found by iRanrirez et al.l (|2012f l 
using optical indices. 

For all stars in our sample we also run the IRFM to 
derive effective temperatures uncorrelated to the spectro- 
scopic analysis. In fact, the IRFM is essentially model 
independent and very little affected by the metallicity 
and surface gravity of each star, the most relevant in- 
gredient being the absolute calibration of the photomet- 
ric systems adopted. Since all these stars are nearby, 
reddening is also not a concern as con firmed by using 
intrin sic Stromgren color calibra tions ()Melendez et al.l 
I2010D . In lCasagrande et al.l (|2010f ) the zero point of the 
IRFM scale was calibrated using solar-twins having Ty- 
cho2 and 2MASS photometry. Now, the availability of 
UBV(RI)c magnitudes allows to check whether this is 
also the case when using instead the Johnson-Cousins 
systeir0. Replacing the Tycho2 with the Johnson- 
Cousins system in the IRFM returns effective temper- 
atures cooler by 26 ± 4 K (cr = 41 K) which is within 
the zero-point uncertainty of our effective temperature 
scale. The same conclusion is obtained restricting the 
analysis to solar-twins onlj0 for which we obtain a me- 
dian/mean T eS of 5787/5786 K and 5762/5750 K using 
Tycho2 and Johnson-Cousins photometry, respectively. 
These differences are mirrored in the color indices of Ta- 
ble |31 those inferred using Johnson-Cousins in the IRFM 
are in fact slightly bluer than obtained via LDRs, by an 
amount which would be almost perfectly offset should 
the effective temperatures increase by 20 K. As expected, 
(V — J)©, (V — H)q and (V — K s )q derived here agree 
extrem ely well with those obtained bv lCasagrande et"al] 
(2010) reverting color— [Fe/H] — T c ff polynomials defined 
over a much wider parameter spaceG3- 

Effective temperatures determined using Tycho2 pho- 
tometry in the IRFM return color indices which agree 

8 For an implementation of the IRFM using WISE magnitudes 
instead of 2MASS see Appendix [X] Here we prefer to use the 
2MASS system because of the better quality data. 

9 I.e. stars having spectroscopic stellar parameters within 1.4<r 
from the solar ones , in accordance with the criterium used by 
IRaimrez efal] 120T3) . 

10 Incidentally, using optical and infrared solar colors from LDRs 
in the aforementioned polynomials returns an average T G ff = 5755± 
22 K. 



almost perfectly with LDRs. These results confirm the 
overall good consistency obtained implementing different 
photometric systems in the IRFM, with systematic un- 
certainties at the level of about 20 K, i.e. ~ 0.01 mag 
in colors. Errors in Table [3] take this into account, by 
adding such zero-point systematic to the uncertainties 
derived analytically from the fits. 

4. THE WISE THROUGH EXCESS OF INFRARED IS MADE 

A FOOL? 

As discussed in Section ® stars having exceedingly red 
indices in WISE (Table [4j were not used to derive so- 
lar colors, although all other photometric quality con- 
straints were satisfied (apart from a few cases having 
in a given band ccf / or signal-to-noise lower than 
10, bands which were always excluded from the analy- 
sis). The unusual colors of these stars are clearly visi- 
ble in the bottom right panel of Figure [TJ The bands 
most strikingly affected are W2 and partly W4, while 
Wl and W3 seem only slightly offset to the red with re- 
spect to the main locus defined by the full sample. This 
sort of signature (Wl — W2 > 0.3) would not be en- 
tirely unexpected if looking at brown dwarfs. In fact, 
WISE's two shortest bands are designed to optimize sen- 
sitivity to this class of objects, by probing their deep 
CH4 absorption band at ~ 3.3/xm (Wl) and the region 
relatively free of opacity at ~ 4.6/im (W2 ) where their 
Planck function approximately peaks (e.g . Wrigh tTet al.l 
[2rJlTl [MahTzer et alj|201ll iKirkpatrick et al.ll20lil) ~ 

To quantify the amount of contamination expected 
from a potential low mass star companion we combine 
a synthetic MARCS (|Gustafsson et all 120081 ) solar spec- 
trum with that of a late type M dwarf (T c ff = 2500 K) of 
the same metallicity, an d assuming a seconda ry to pri- 
mary radius ratio of 0.2 (|Bovaiian et al.|[2012b| ). we con- 
clude that the effect on the color indices shown in Figure 
Q] would be of order 0.01 mag. Therefore, even if the 
spectral features of a brown dwarf could account for the 
red Wl — W2 index we observe, the overwhelming flux of 
the primary makes this solution not viable even for an M 
dwarf, in addition to the fact that it would still be diffi- 
cult not to affect W3 . Neither the alignment/confusion 
with extragalactic sources (which would be considerably 
fainter than our objects, see below) or cool (sub)-stellar 
objects is likely The angular resolution of WISE passes 
from 6.1 arcsec in Wl to 12 arcsec in W4; using the 
higher resolution of 2MASS, none of the targets discussed 
here has more than one counterpart within 12 arcsec. 
All anomalous sources have Wl brighter than 7.6 mag 
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(5.9 mag if considering the reddest V — W2 > 1.8). Us- 
ing this constraint with the previous synthetic model 
(T G ff = 2500 K) would imply the presence of cool M 
dwarfs closer than 7.5 (3.5) pc and brighter than V ~ 17 
(15.5) mag. Using ins tead a synthetic brow n dwarf spec- 
trum (T e fj = 1000 K, iBurrows et alJl2~006l ) and adopting 
R = 0.1 Rq the above estimates would change into dis- 
tances closer than 1 (0.5) pc and V magnitudes brighter 
than ~ 21 (19.5) mag, thus making extremely unlikely 
the superposition of brown dwarfs to our solar-like stars. 

Interpreting the color anomaly as mid-infrared emis- 
sion could hint toward the presence of warm/hot debris 
disks, even though this class of object is thought to be 
very rare compared to most of the known cold Kuiper- 
bclt type disks, especially around old ma in sequence 
stars (e.g. iBrvden et all [2001 lWvattll2008l ). In Figure 
Q] it is interesting to include for comparison HD 69830 
(HIP 40693), a - 2 - 4 Gyr old, solar metallicity 
dwarf known to host a warm disk closer than ~ 1 AU 
(|Beichman et al.l I2005D as well as three N eptune-mass 
planets within 0.6 AU (jLovis et al.l I2006D . Although 
HD 69830 has a T e g somew hat cooler than the bulk of 
our stars (|Sousa et al.l 12008). it shows a clear excess in 
W2 and 1U4, while Wl and W3 are barely affected. 
A somewhat similar trend is also shown by HD 72905 
(HIP 42438), a ~ 0.2 Gyr old, sola r- like star surrounded 
by hot dust (jBeichman et al.l l2006;FL 

However, producing emission in W2 while keeping the 
other two contiguous bands essentially unaffected re- 
quires major fine tunings, thus rendering also the disk 
interpretation difficult. At first, interpreting all anoma- 
lous WISE colors as suspected and/or poor photometry 
seems difficult because of the various quality constraints 
imposed (Section [5J , although a number of W2 excesses 
have a saturated pixel fraction higher than usual (c.f. Sec- 
tion [2l and more st r ingen t cuts do alleviate the problem. 
Kennedy fc Wv att (2012|) have conducted a throughout 
study of stars in the Kepler field using WISE photom- 
etry to identify disk candidates, and concluded in fact 
that spurious detections are less likely when using pho- 
tometric measurements with source extraction x 2 smaller 
than 2, 1.5, 1.2, 1.2 in Wl — A, respectively. The require- 
ment ext = we adopt is in fact less stringent (Section 
[5]); the tighter % 2 constraints listed above are satisfied 
by 98% (Wl), 80% (W2), 53% (W3), 93% (WA), of the 
stars used in Section H and by 100% (Wl), 50% (W2), 
50% (W3), 95% (WA) of the stars in Table H This sug- 
gests that any excess we see in WA might be real, but 
it casts some doubts on W2. Adopting these constraints 
the number of outliers in Table 0] considerably reduces, 
although a number of them still remains: HIP 109110, 
objects with excess in W2 only (HIP 38228, HIP 88194 
HIP 100963) and WA only (HIP 38072 and HIP 66885). 

For these six stars (plus HD 69830 and HD 72905), 
photospheric models tailored at the measured spectro- 
scopic parameters are compared to fluxes derived from 

11 HD 69830 satisfies all photometric quality requirements listed 
in Section [2] apart from missing 2MASS counterpart within 3 arc- 
sec. Notice though that this star has proper motion of ~ 1 arc- 
sec/year. HD 72905 was in our initial sample of Section [2] but 
discarded because ext = 1. This flag implies that the profile fit of 
the photometry is not optimal (namely x 2 = 3.6 in W3), although 
it is still not associated with any 2MASS Extended Source and all 
other photometric quality requirements of Section [2] are satisfied. 



the adopted photometry (Figure [3]) . The absolute cali- 
bration of those model fluxes is done by forcing them to 
return synthetic V magnitudes that match the observed 
ones. We also tested on the full sample of stars in Section 
[2] that the mean difference in physical fluxes is 0.25% and 
never exceeds 1.5% if doing instead the absolute calibra- 
tion using angular diameters obtained from the IRFM. 
These differences are essentially indistinguishable on the 
scale of the plot, and are taken into account when com- 
puting the flux uncertainty associated to each photomet- 
ric measurement. Magnitudes in the BV(RI)c JHK S 
system are converted into fl uxes using the same ab - 
solute calibration adopted in [Casagran de et "all |2010), 
which has an intrinsic uncertainty of order 1 — 2%. Be- 
cause of the aforementioned difference between using V 
magnitudes or angular diameters, we have increased the 
global flux error to a conservative 3%. Similarly, for 
the WISE system, we have adopted the a bsolute cali- 
bration and errors from lJarrett et all (|2011[ ). further in- 
creasing the latter by 1%, which also takes into account 
the possible zero-point offset discussed in Appendix 
Observed magnitude errors are then added to the afore- 
mentioned uncertainties regarding the absolute flux of 
each band. As already expected from color indices, any 
difference with respect to photospheric models in W2 
and WA is significant, and it does not stem from un- 
certainties on the flux scale, nor in the observed mag- 
nitudes. The advantage of fitting photospheric models 
instead of using color indices is that we are now able to 
better quantify the observed anomalies. The same com- 
parison with synthetic spectra have also been done for all 
other stars in our sample not showing anomalous WISE 
colors, and indeed there are no mismatches between syn- 
thetic spectra and observations, thus validating the over- 
all flux scale we adopt and also excluding major model 
deficiencies in those ban ds for solar-like stars (but see 
iKennedv &: W vatt 2012, for possible model inaccuracies 
at cooler T c s). 

As we already discussed, adding a cool companion does 
not modify the energy distribution in a way which is 
able to explain the observations, apart from HIP 109110, 
which longward of the J band shows fluxes systematically 
higher than predicted. This likely indicates the presence 
of a cool companion (which we are able to fit with a 
model having T c s ~ 4000 K), an inter pretation which is 
consi stent with its suspected binarity (jFrankowski et all 
l2007h and with the linear t rend observed in its radial ve- 
locity (jNidever et al.ll2002h . This infrared excess is then 
further con firm by the Spitzer-IR AC and -MIPS mea- 
surements ([Carpenter et al.|[2008f ) shown in Figure [3] as 
open squares. 

The W2 measurements for HD 69830 and HD 72905 
do not pass the more stringent requirements we impose 
on the source extraction. The spurious nature of W2 
photometry for these two stars is then confirmed by 
the absence of excess in other measu r ements at simi- 
lar w avelengths (|Carpenter et alj 1200 8: Bcichma n et al.l 
|20H . Note tho ugh, for HD 69830 the . IRS, Syitzer-MIPS 
and IRAS daXa. ([Beichman et al.ll2011l ) confirm the excess 
we see in WA. Thus, we are inclined to regard any W2 
excess among our stars as artificial, even if the source ex- 
traction is fine; this is further supported by the fact that 
in Figure [3] the W2 phot ometry of H IP 100963 does not 
agree with Svitzer-IRAC ([Ca,rventer et al.ll2008f) . No ad- 
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ditional measurements around 4.6/im exist for HIP 38228 
and HIP 88194, but from the previous discussion, and be- 
cause these excesses seem rather challenging to interpret 
when contiguos bands agree well with photospheric mod- 
els, we conclude that also their nature is likely spurious. 

Finally, for HIP 38072 and HIP 66885 the deviation 
from photospheric models starts only in WA. Despite 
these two stars being the faintest among those in TablelH 
so far our adopted quality constraints have been enough 
to discard unreliable WA measurements, and what we 
see could be indeed the signature of debris disks around 
these two stars. For HIP 38072 we have a measured 
flux of 10.7 ± 0.7 mJy and a photospheric prediction of 
6.7 ± 0.2 mJy, thus resulting in an excess ratio of 1.6 
with a 5cr significance, while for HIP 66885 the measured 
flux is 8.7 ± 0.8 mJy versus a photospheric prediction of 
7.0 ± 0.2 mJy, the excess ratio being 1.2 at 2a. Using 
the absence of emission in W3 to constrain their tem- 
perature, we are able to easily fit these excesses with 
a black-body, but measurements at other wavelengths 
are clearly required to confirm or discard the presence 
of any disk. Should these two detections be confirmed, 
and using the variance of the the binomial distribution to 
derive a realistic error ba r for such a low number statis- 
tic (e.g. lBevingto"nlll969D the occurrence rate of debris 
disk at 22/Ltm from our solar-like sample would thus be 
3 ± 2%, in good agr eement with the ~ 4% estimated by 
iTrilling et atl ()2008h at 24/im. 

5. CONCLUSIONS 

The uncertainty on the zero point of the effective tem- 
perature scale has been addressed using one of the most 
stringent photometric contraints available, i.e. the in- 
frared colors of Sun, here determined in a model indepen- 
dent way via LDRs. Such analysis leads to an excellent 
agreement with the indic es derived using the T e g scale of 
ICasagrande et al.l (|2010D . thus confirming its accuracy. 
Effective temperatures have also been determined imple- 
menting the IRFM in the WISE system, thus validating 
the overall consistency between the 2MASS and WISE 
absolute calibration. This work also shows the impor- 
tance of using solar twins for the absolute calibration 
of photometric quantities, something which is getting in- 
creasingly more important in the era of large all-sky pho- 
tometric surveys. 

However, while 2MASS magnitudes are very well 
suited for the purpose of determining T c g -once stars 
affected by binarity are removed and/or using the full 
quality and flag information available in 2MASS to dis- 
card bad photometry- this does not necessarily hold for 
the WISE data. Special attention on WISE's photomet- 
ric quality flags and source extraction information must 
be paid, yet a number of excess emissions in W2 seem ar- 
tificial. In this respect 2MASS magnitudes lie in a "sweet 
spot" , enough in the red to sample the Planck tail, yet 
largely unaffected by contamination and/or flux excess, 
that being real or spurious. 

At this stage it is still unclear if all WISE mid-infrared 
anomalies found are stemming from bad measurements 
or they are rather associate to real physical phenomena. 
Either cases being possible, those stars are clearly not 
representative of the Sun and have not been used to 
derive its colors. It is extremely difficult to interpret 
in a consistent manner objects showing intense excess 



Table 4 

Stars having V - W2 > 1.65 or V - WA > 1.80. 

Bands that do not satisfy all of the quality 
constraints discussed in Section [2] are written in 
italics. The x 2 of the source extraction in each 
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in W2 only (and in fact, comparison with independent 
measurements confirms the spurious nature), while on 
the contrary it seems genuine for WA. This latter ex- 
cess could be the signature of warm/hot debris disks, 
the best candidates from our sample being HIP 38072 
and HIP 66885. Data at longer wavelengths are clearly 
needed: should HIP 38072 be confirmed, it would be the 
first solar-analog/twin (T S g CC = 5839 K, \ogg = 4.53 dex, 
[Fc/H] = 0.06 dex) found to host such a debris disk. This 
star is also relatively young (2.4 Gyr) and has ~ 4 times 
more lithium than the Sun dBaumann et^[20T0t) . and 
it is included in our HARPS radial velocity monitoring, 
thus making it a potentially interesting target to gauge 
new insights into the planet-disk interaction. 

Note added in proof. The reason for the anomalously 
bright 1^2 magnitudes likely stems from problems in the 
profile fitting photometry for bright sat urated stars, see 
http : / / wise2 .ipac.caltech.edu/docs/release/ allsky / expsup / 
sec6_3c.html 

We thank an anonymous referee for a prompt report 
and G. Kennedy for pointing out the WISE bias in the 
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Figure 3. Observed fluxes vs. photosphcric models (from Castclli & Kurucz 2004) for stars still having anomalous WISE colors after 
imposing more stringent requirements on the source extraction. Fluxes — F v (3y) — are derived from optical (cyan), 2MASS (green) and 
WISE (red) magnitudes as described in text: filled (open) symbols identify WISE photometric measurements that do (not) satisfy the 
requirements discusse d in Section fj] includ i ng the stricter \ 2 extraction threshold. O pen squares are independent flux measurements 
available in literature (Carpenter ct al. 2008; Plavchan ct al. 2009; Bcichman ct al. 2011). 
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Figure 4. AT c ff (2MASS — WISE) when implementing one or another system in the IRFM. Filled squares are stars having V — W2 and 
V — WA redder than the thresholds discussed in Section [2] and are not used to derive the mean difference (continuous line). Left panel: 
using in the IRFM at least two bands among W1W2W3W4. Right panel: when restricting to at least two bands among WIW2W3. 
Gray shaded area is the systematic offset (±36 K) allowed by the uncertainty in the WISE's absolute calibration. 
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APPENDIX 

THE INFRARED FLUX METHOD USING WISE PHOTOMETRY 

The availability of 2MASS and WISE photometry for most of the stars in our sample permits to run the IRFM using 
these two systems separately, and to check th at consistent effect i ve temperat ures are derived. The implementation 
of the method is identical to that described in [Casagran de et al.l (|2006l l20Tol ). i.e. for each star the bolometric flux 
is recovered using all available broad-band optical/near- infrared/mid- infrared colors, while infrared fluxes needed to 
derive T e g are now computed using 2MASS or WISE photometry, respectively. 

For solar-type stars, the main driver in setting the zero-point of the T B ff scale is the a bsolute calibration of infrared 
bands. We have already discussed and tested that of 2MASS (|Casagrande et al1l2010L and references therein), thus 
focusing on the WISE system her e. We adopt t he W l W2 W3 WA relative system response curves H and physical 
monochromatic fluxes -F^iso) from lJarrett et all (|2011f ). The latter are built on the same absolute b asis established for 
the Spitzer Space Telescope, ultimately constructed on the "Cohen- Walker- Witterborn " framework (jCohen et al.ll2003l . 
and ref erences therein) a nd tied directly to the absolute mid-infrared calibrations by the Midcourse Space Experiment 
(MSX, P rice et aLll2004f) . These WISE fluxes have an expected overall systematic uncertainty of order ~ 1.5% and 
define the Vega zero-magnitude attributes upon which the effective temperatures we derive in this system directly 
depend. For each star we used the photometric constraints discussed in Section [2] and derived T c g if at least two WISE 
bands were simultaneously available. WA magnitudes have the largest photometric errors (Section [2]) thus showing 
the weakest correlation with T ff ; indeed the scatter in the comparison with the effective temperatures obtained using 
2MASS magnitudes reduces when WA is not implemented in the IRFM (Figure U). 

The maximum zero-point uncertainty in the effective temperatures derived using WISE photometry can be easily 
estima ted by increasing or decreasing the adopted absolute calibration according to the error reported in lJarrett et al.l 
(|2011l ). As expected, for the stellar parameter space covered here, the effect is essentially a constant offset of ±36±0.4 K 
a = 3 K). Notice that our reference 2MASS effective temperatures also have a zero-point uncertainty of order 20 K 
Casagrande et al.ll20l6T) . The 25 - 30 K mean difference when using 2MASS or WISE magnitudes in the IRFM is thus 
fully consistent within the uncertainties, and it would disappear if the WISE absolute calibration had been decreased 
by about 1% (or conversely, the 2MASS absolute calibration increased by the same amount). Such a nice agreement is 
not entirely unexpected, since the 2MASS and WISE absolute calibration are built within the same " Cohen- Walker- 
Witterborn " framework: nevertheless, the fact that the 2MASS absolute calibration has been independently verified 
using solar-twins confirms that the accuracy of the infrared absolute calibration extend also to the mid-infrared regime 
probed by WISE. 



http: / / wise2.ipac.caltech.edu/docs /release / allsky / expsup /sec4_4h.html 
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Table 1 

(V — K s )q colors inferred from LDR Measurements. AT* is the number of 
stars used for the fit, (jfi t is the standard deviation of the fit minus data 
residual, and cr BB is the standard deviation of the color inferred from the nine 
reflected Sun-ligh asteroid observations used for solar reference. 
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